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ABSTRACT 

The discovery of rings around extrasolar planets (“exorings”) is one of the next breakthroughs in 
exoplanetary research. Previous studies have explored the feasibility of detecting exorings with present 
and future photometric sensitivities by seeking anomalous deviations in the residuals of a standard 
transit light curve fit, at the level of ~ 100 ppm for Kronian rings. In this work, we explore two much 
larger observational consequences of exorings: (1) the significant increase in transit depth that may 
lead to the misclassification of ringed planetary candidates as false-positives and/or the underestima¬ 
tion of planetary density; and (2) the so-called “photo-ring” effect, a new asterodensity profiling effect, 
revealed by a comparison of the light curve derived stellar density to that measured with independent 
methods (e.g., asteroseismology). While these methods do not provide an unambiguous detection of 
exorings, we show that the large amplitude of these effects combined with their relatively simple ana¬ 
lytic description, makes them highly suited to large-scale surveys to identify candidate ringed planets 
worthy of more detailed investigation. Moreover, these methods lend themselves to ensemble analyses 
seeking to uncover evidence of a population of ringed planets. We describe the method in detail, 
develop the basic underlying formalism and test it in the parameter space of rings and transit config¬ 
uration. We discuss the prospects of using this method for the first systematic search of exoplanetary 
rings in the Kepler database and provide a basic computational code for implementing it. 

Subject headings: Techniques: photometric Occultations — Methods: analytical — Planets and 
satellites: rings 


1. INTRODUCTION 

Since the discovery of the first transiting planet (Char- 
bonneau et al. 2000; Henry et al. 2000), planetary eclipses 
have emerged as powerful tools for characterizing exo¬ 
planets. Numerous novel methods have been devised us¬ 
ing transits to identify non-conventional planetary prop¬ 
erties, such as oblateness (Carter & Winn 2010; Leconte 
et al. 2011), magnetic bow shocks (Vidotto et al. 2010), 
exomoons (Kipping 2009; Heller et al. 2014), and exo¬ 
planetary rings or “exorings” (Barnes & Fortney 2004; 
Ohta et al. 2009; Tusnski & Valio 2011; Mamajek et al. 
2012; Kenworthy & Mamajek 2015). 

With the advent of new instruments surveying the sky 
for transits with precise photometry (PLATO, Rauer & 
Cat ala 2011; EELT, Guyon et al. 2012; GMT, Johns 
et al. 2012; TESS, Ricker et al. 2014 and JWST, Beich- 
man et al. 2014), there is great potential for discoverying 
new and unconventional exoplanetary phenomena in the 
coming decade. 

The discovery of exorings would be particularly in¬ 
teresting. All of the solar system’s giant planets have 
rings and at least one, Saturn, has sufficiently extended 
ring systems that are sufficiently extended to produce ob¬ 
servable signatures with current/future instrumentation 
(Barnes & Fortney 2004; Ohta et al. 2009). The discovery 
and characterizarion of exorings could shed light on im¬ 
portant planetary processes, such as planetary and moon 
formation (Mamajek et al. 2012) and planetary interior 
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structure (Schlichting & Chang 2011). 

Barnes & Fortney (2004) developed the first theoreti¬ 
cal model of exoring transits. They showed that, pro¬ 
vided photometric sensitivities and time resolution of 
1 —3xl0 -4 and ~ 15 minutes, light curve residual analy¬ 
sis could be used to resolve and characterize rings. More 
recently, Ohta et al. (2009) showed that the presence of 
rings could produce spectroscopic signatures, detectable 
by ground-based telescopes, for sensitivities < 0.1 m 
s _1 . Independently, Tusnski & Valio (2011) have tested 
a model of ring transits showing that the properties of 
hypothetical rings can be reliably recovered. 

Despite these theoretical advances, there have been no 
systematic surveys for exorings with archival photome¬ 
try. One reason for this is the assumption that a practi¬ 
cal ring survey would employ the same technique as that 
for their discovery and characterization, namely, the de¬ 
tailed analysis of light curves via fits to complex ring 
transit models. 

Here, we present a novel method for perfoming system¬ 
atic searches for exoring candidates. Our method relies 
only on the measurement of the basic transit parame¬ 
ters (i.e. depth, duration), avoiding the need for com¬ 
putationally expensive fits of a transit ring model. The 
method exploits simple analytical formulae and hence is 
well-suited for performing searches in large photometric 
databases. 

2. RING TRANSIT GEOMETRY 

One of the most powerful properties of planetary tran¬ 
sits is the wealth of information they can provide by sim¬ 
ply monitoring variations in stellar brightness. 

Assuming that the transiting object is spherical and 
other conventional conditions, Seager & Mallen-Ornelas 
(2003) showed that the planetary radius p = 


2 



Fig. 1.— Schematic ring and planetary transit geometry and light curve. Sizes are not at scale. Transit depth (S) and total duration 
(T 14 ) are larger when an exoring is present, whereas the duration of full transit (T 23 ) is lower. 


scaled orbital semimajor axis, a/i?*, and impact param¬ 
eter, b, can be derived from three basic observables: 


Ar p = nR 2 p + 7r[r 2 (/ e ) - r 2 (/i)] R 2 , (2) 


1. transit depth, 5 = ( F a — F)/F 0 (where F and F 0 
represent the in-transit and out-of-transit stellar 
fluxes, respectively); 

2. first-to-fourth contact transit duration, Tj 4 ; 

3. second-to-third contact transit duration, T 23 . 

Figure 1 schematically depicts the definition of these 
quantities and the significant differences imposed by the 
presence of planetary rings. 

Our basic ring transit model relies on three basic as¬ 
sumptions. 


1. The planet and the star are spherical. 

2. Rings are uniform and scatter/absorb light only be¬ 
tween radii Ri = fiR v and R e = f e R P (constant 
normal optical depth r). 

3. Diffractive forward scattering in ring particles does 
not modify the basic transit parameters (Barnes & 
Fortney 2004). 


Under these assumptions, the transit depth of a ringed 
planet is given by the ratio: 


5 = 


Ar p 

A* 


Ar p 

ttRI 1 


( 1 ) 


where Ar p and A* are the effective projected ring-planet 
and stellar area. 

The ring area is computed with the analytical expres¬ 
sion (see Appendix A for details) 


where r is the effective ring radius and is given by 


i f 2 COS 2R —1 /cOSlR>l 

f 2 cos zr |-arcsin(y) + 

^ otherwise 

- — arcsin(y/cos ir) 

( 3 ) 

Here, zr is the projected ring inclination (zr = 90° if 
the ring is edge on) and y = ^ f 2 — 1/(/suizr) is an 
auxiliary variable. The term /?(coszr) = 1 — e -T / coslR 
accounts for the ring’s effective absortion of stellar light 
(Barnes & Fortney 2004). 

To calculate ring transit durations, T 14 and T 23 , we 
need to compute the planet center horizontal coordinate, 
x, where the the external ring or the planetary disk in¬ 
tersects the stellar limb at a single point (see Figure 1). 
Four x values, the contact positions, x% (z = 1,2, 3,4), 
fulfill this condition. 

In the case of a non-ringed spherical planet, contact 
positions are given by 


< 4 > 

where p = R p /R+ and the “+” sign applies to contacts 1 
and 4 and the sign to contacts 2 and 3. 

If the contact point is on the edge of the external 
ring, then this point satisfies a set of non-trivial alge¬ 
braic/trigonometric equations, whose solution is not ex- 
pressable in a closed form (see Appendix B). Combining 
several approximations, we have found that in a wide 
range of ring and transit configurations, the following 
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analytical formula provides the value of ( Xi/R +) with a 
relative uncertainty no larger than a few percent: 


± Acos 6 »^ « 1 - A 2 (sin0 =f b/A) 2 (l - B 2 /A). 

( 5 ) 

Here, A = f e p and B = f e p cos iu, are the ring pro¬ 
jected semimajor and semiminor axes, and the signs now 
correspond to temporal ordering. The upper signs corre¬ 
spond to the “leading” contacts (contacts 1 and 3) and 
the lower ones to the “trailing” contacts (contacts 2 and 

4 ). 

The positive solution to both, Eqs. (4) and (5), corre¬ 
sponds to contacts 3 and 4 and the negative solution to 
contacts 1 and 2 . 

Once the contact positions are calculated, the total 
duration of the transit X 14 and the duration of full transit 
T 23 may be estimated using (Sackett 1999), 


a .. . 

— sin 1 sin 
rik 

a . . . 

—— sin 1 sin 
R-k 



X4 

Xi 


R* 

X 3 

X 2 

R* 

R* 


( 6 ) 

( 7 ) 


where i and P are the inclination and period of the plan¬ 
etary orbit. 

3. RING EFFECT ON OBSERVED PLANETARY RADIUS 

Ring transits could have a significant effect on the tran¬ 
sit depth, S. To first order, <5 provides the value of the 
apparent or observed planetary radius, p 0 bs: 

Pob s = VS. ( 8 ) 

Thus, a significant overestimation of 6 will also produce 
a substantial overestimation of the observed radius. In 
the absence of any other evidence for the presence of 
exorings, an anomalously deep transit may lead to the 
misclassification of a ringed planet as a false-positive or 
yield a gross underestimation of the planetary density. 

To illustrate this, hereafter we assume the reference 
case of a planet with the same radius as Saturn R p = 
O.O836i?0 and rings of similar size and properties, i.e., 
/i = 1.58 (inner edge of the A-ring) and f e = 2.35 (outer 
edge of the B-ring). We assume for simplicity a value of 
r = 1 (the opacity of B-ring ranges from 0.4-2.5, Murray 
& Dermott 1999). We also assume that the planet orbits 
a solar-mass star in a circular orbit at a = 1 AU. The 
generalization of these results to other planetary, ring, 
and orbital parameters is trivial. 

In Figure 2, we plot the ratio of observed to true plane¬ 
tary radius as obtained with Eqs. (1-3) for our reference 
planet. 

In the extreme case of face-on rings, Eqs. (1-3) simplify 
to give 


lir Sn 0 P° hs = P\J l + ( l - e T )(fe-f ?)• ( 9 ) 

jr->90° V 

The transits of a Saturn-like ringed planet are up to 
~3 times deeper than these expected for a spherical non- 
ringed one. These deep transits will be interpreted as 




Fig. 2.— Upper panel: ratio of observed to true planetary radius 
as a function of ring projected inclination (cos i), assuming differ¬ 
ent values for the ring normal opacity (r). Lower panel: probability 
distribution of this ratio assuming completely random ring orien¬ 
tations (blue histogram) and a more concentrated distribution of 
planetary obliquities (red dashed line). For illustration purposes, 
we show, in the inset diagram, the distribution of random axis ori¬ 
entations used to compute both distributions. Each point represent 
the positions of the north and south pole of the planet for a given 
random orientation. 

being produced by a planet ~1.7 times larger. Addi¬ 
tionally, if independent estimations of its mass were also 
available, then the density of the planet would be under¬ 
estimated by a factor of ~5. Thus, instead of measuring 
Saturn’s density ~0.7 g cm -3 , this planet would seem 
to have an anomalously low density of ~0.14 g cm -3 . 
Even under more realistic orientations (coszr ~ 0.2), the 
observed radius will be ~ 20 % larger and the estimated 
density almost a half of the real one. 

To assess the effect of rings on the observed planetary 
radius for an ensemble population of planets, we have cal¬ 
culated the probability distribution of the ratio (p Q bs/p) 
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assuming a uniform random ring orientation (obliquity 
and azimuthal angle, as defined in Ohta et al. 2009). 
The result is plotted in the lower panel of Figure 2 (blue 
histogram). The probability distribution corresponding 
to predominantly low obliquities (following a Fisher dis¬ 
tribution in the case of a concentration parameter k ~ 8 
and 2-cr obliquity dispersion of ~ 30°) is also shown for 
comparison (red dashed line). 

In the case of uniform random obliquities, more than 
50% of the orientations lead to overestimations in plane¬ 
tary radius greater than ~50%. This implies that with no 
other clues for the existence of rings around those plan¬ 
ets, the bulk density of more than half of them would be 
underestimated by up to a factor of 3. With a more con¬ 
centrated distribution of obliquities, the resulting radius 
distribution peaks at a null-effect of ( p 0 bs/p ) ~ 1 but 
still exhibits considerable dispersion, with >50% of the 
cases having observed radius anomalies > 30% (density 
overestimated by a factor > 2 ). 

Deep transits are used as a criterion for flagging poten¬ 
tial false-positives in photometric surveys (Batalha 2014; 
Burke et al. 2014). If a large unseen population of ringed 
planets exists, then we could have many potential candi¬ 
dates buried among misclassified false-positives. There¬ 
fore, we suggest that revisiting false-positive transits re¬ 
jected with this criterion could potentially lead to the 
discovery of the first exoring population. 

4. THE PHOTO-RING EFFECT (PR-EFFECT) 

Under the simplified assumptions of the model pre¬ 
sented in Section 2, an analytic solution for the orbital 
semimajor axis and impact parameter can be obtained 
from the basic transit parameters (Seager & Mallen- 
Ornelas 2003): 



obs 


^obs 


P (5 1 / 4 
2^ (T? 4 - r 2 2 3 ) 1/2 


1/2 


rr 2 _ rp2 
1 14 1 23 


( 10 ) 

( 11 ) 


where we also assume Ti 4 ,T 23 -C P, which is suitable 
for planets with a i?*, i.e., those worlds with improved 
chances for exorings. 

Kepler’s third law relates the orbital period, semimajor 
axis, and stellar mass, such that we can calculate the 
mean stellar density: 


P*,obs 


3tt (a/-R*)g bs 
G P 2 


( 12 ) 


In the case of a transiting spherical planet, the ob¬ 
served value of (a/i?*)obs, and hence p*, 0 bs, are accurate 
estimates of their true values. If the planet, however, has 
a ring, then S , T 14 and X 23 will not be related by Equa¬ 
tions ( 10 ) and ( 11 ), and the observed quantities will differ 
from the true ones. 

A comparison of the observed density p* j0 bs with an 
independent measurement, p* i true, such as that provided 
by stellar models, asteroseismology (Huber et al. 2013), 
or transits of other planetary companions, could indicate 
the existence of exorings. 

In the upper panel of Figure 3, we show contours in the 
plane of the projected ring orientation of this so-called 


Photo-Ring Effect 
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Fig. 3.— Upper panel: contour plot of the PR effect. The crosses 
correspond to actual values of the projected inclinations for a sub¬ 
sample drawn from a low obliquity population (red dots in the inset 
plot of the lower panel). Lower panel: probability distribution of 
the PR effect for completely random orientations (blue) and for a 
more concentrated distribution of obliquities (red). 


PR-effect. 

Since T 14 (T 23 ) is larger (smaller) when a ring is present 
(see Figure 1), the observed semimajor axis a Q bs could 
be over- or underestimated. For large tilts (6 ~ 90°) and 
inclinations (coszr ~ 0 ), the transit duration will be the 
same as that of a non-ringed planet. Consequently, the 
increased transit depth 5 will be the dominant effect in 
Eq. (10) and the observed density will be overestimated 
(upper left region in Figure 3). 

For most values of the ring’s projected inclination and 
tilt, however, the transit duration will be modified to 
a larger extent than the depth, causing the semimajor 
axis and density to be underestimated. Therefore, rings 
tend to produce a negative PR effect (underestimation of 
stellar density) and this could be used to distinguish it 
from other asterodensity profiling (AP) effects discussed 
in Kipping (2014). 
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We can use Eqs. (4)-(7) to calculate an analytical ex¬ 
pression for the maximum PR effect expected for a given 
external disk radius and impact parameter: 

«/- 3/2 (l-6 2 )- 3/4 (13) 

\P*,true / max 

For our reference case, f e =2.35, 6 = 0, and the max¬ 
imum value of the PR effect is log 10 (/9* iO bs/p*,true) ~ 
—0.6. This corresponds to an underestimation of the 
stellar density by a noticeable factor of ~ 4. 

We have calculated the probability distribution of the 
PR effect assuming similar priors to those used earlier. 
The results are shown in the lower panel of Figure 3. 

In the uniform case, the PR distribution is strongly 
peaked around a negative value, log 10 (p* i0 bs/p*,true) ~ 
—0.2. In contrast, if we consider obliquities no larger 
than ~30°, then the peak shifts to log(p* i0 bs/p*,true) ~ 
—0.45, corresponding to a notable difference by a factor 
of ~ 3 between the observed and true density. 

The only other AP effects which can cause compara¬ 
bly large deviations are the photo-eccentric (PE) and the 
photo-blend effects, for which the former tends to cause a 
positive AP deviation and the latter a negative (Kipping 
2014). In the case where blended companions can be ex¬ 
cluded (e.g., through high-resolution imaging), then the 
ensemble distribution of the PR-effect should therefore 
be distinguishable from the other AP effects. 

5. DISCUSSION AND CONCLUSIONS 

In this Letter, we have presented a novel method for 
identifying exorings. Our method does not require a com¬ 
plex fit of transit light curves, instead relying on simple, 
analytic, and computationally efficient numerical proce¬ 
dures. 

Our technique exploits the substantial impact that 
rings produce on the transit depth and duration, effects 
that can be independently confirmed by the measure¬ 
ment or estimation of other relevant astrophysical prop¬ 
erties (e.g., stellar and planetary density). 

Interestingly, the two effects we seek (anomalous tran¬ 
sit depths and PR-effect) are complementary with re¬ 
spect to the orientation of the ring plane. For large in¬ 
clinations and obliquities (face-on rings), the effect on 
transit depths is significant whilst the photo-ring effect 
is negligible. Alternatively, if rings have relatively low 
obliquities (egde-on rings), then the PR-effect will be 
considerable but the depth anomaly small. 

Accordingly, three basic complementary strategies are 
proposed to identify exorings among the confirmed tran¬ 
siting exoplanets and candidates. 

1. Search the confirmed transiting planets with 
anomalously low densities. 

2. Search transiting objects that have been tagged 
as false positives due to anomalously large transit 
depths. 

3. Search transit signals for which a negative AP effect 
is observed. 

Besides aiding observers seeking interesting individual 
systems, the effects described here are well-suited for in¬ 
ferring a population of ringed planets. Ensemble studies 


using AP are in their infancy, but Sliski & Kipping (2014) 
recently conducted the first such analysis on a sample 
of 41 single Kepler planetary candidates with asteroseis- 
mically constrained host stars. In this limited sample, 
largely focused on short-period planets, the authors con¬ 
clude that the 31 dwarf stars yield a broad AP distri¬ 
bution about zero, consistent with the expected photo¬ 
eccentric variations, whereas as the objects associated 
with giant stars are likely orbiting different stars alto¬ 
gether. A larger sample, including longer-period planets, 
would provide the opportunity to seek the expected offset 
due to the PR-effect too. This would require a reliable 
independent measure of the stellar density for stars too 
faint for asteroseismology, perhaps using methods such 
as “flicker” (Kipping et al. 2014a). 

Low-density planets are also interesting targets when 
seeking exoplanetary rings. Over 10% of the already 
confirmed planets have estimated densities 5 . Most of 
them are planets larger than Neptune, and a non- 
negligible fraction have anomalously low densities below 
~0.3 g cm -3 . Although several successful explanations 
have been devised for reconciling observed low densities 
with planetary interior and thermal evolutionary mod¬ 
els (Miller et al. 2009), other anomalies still remain and 
may be worth a further analysis along the lines suggested 
here. 

With the exception of Kepler- 421b (Kipping et al. 
2014b), all of the confirmed planets and most of Kepler 
candidates are inside the so-called snow line. Although 
icy rings, such as those observed around Saturn and So¬ 
lar System giant planets, seem unlikely inside this limit, 
the existence of “warm” rocky rings at distances as short 
as ~0.1 AU, are not dynamically excluded (Schlichting 
& Chang 2011). 

We stress that the method presented here is comple¬ 
mentary to the methods developed to discover exorings 
through detailed light curve modeling (Barnes & Fort¬ 
ney 2004; Ohta et al. 2009; Tusnski & Valio 2011). As 
explained earlier, the role of these methods will be very 
important once a suitable list of potential exoring can¬ 
didates is created. It is, however, also important to note 
the great value of light curve models developed under 
the guiding principle of computational efficiency (semi- 
analytical formulae, efficient numerical procedures, etc), 
such as the basic models presented here. 

Our simple technique is suitable for surveying en¬ 
tire catalogs of transiting planet candidates for exor¬ 
ing candidates, providing a subset of objects worthy 
of more detailed light curve analysis. Moreover, the 
technique is highly suited for uncovering evidence of a 
population of ringed planets by comparing the radius 
anomaly and PR-effects in ensemble studies. To aid 
the community, we provide the publicly available code at 
http://github.org/facom/exorings to simulate the novel 
effects described in this work. 
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Fig. 4.— Geometrical constructions used to compute ring area (left panel) and contact position (right panel). For contact position, a 
rotated reference system where the rings are horizontal considerably eases the analytic derivations. 
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APPENDIX 
AREA OF THE RING 

The area obscured only by the ring is computed substracting the areas of the ellipse and circle sectors limited by 
interesection points Ei and E 2 (Figure 4): 


Aring — 2(A e — Ac). (Al) 

The area of the circle and ellipse sectors follow from the Cavallieri’s Principle 6 : 

A c = i AO R 2 (A2) 

A e = E AB, 

where A and B are the apparent semimajor and semiminor axes of the ring, AO is the angle subtended by the sector, 
and AE is the eccentric anomaly difference between Ei and E 2 . 

The length of the segment joining E]E 2 is 7 : 



A 2 - R? 
A 2 - B 2 


(A3) 


Simple trigonometrical relationships and the parametric equation of the ellipse provide us with the expressions for 
AO and AE: 


AO = 2 arcsin 
AE = 2 arcsin 



Inserting Equations (A5) into A3, the areas are finally given by 


(A4) 


6 A. Bogomolny, Interactive Math¬ 
ematics Miscellany and Puzzles 

http: //www.cut-the-knot.org/Generalization/Cavalieri2.shtml, 


Accessed 2015 January 21 

7 Weisstein, Eric W. MathWorld. 

http: //mathworld. wolfram.com/Circle-Ellipselntersection.html. , 
Accessed 2015 January 21 
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A c = R arcsin cos *r 


A e = f 2 cos *r R 2 arcsin 




f sin *r I 

vz!zr 

/ sin * R 


From these areas and Eq. (Al), the effective ring radius in Eq. (3) follows trivially. 

RING CONTACT POSITIONS 

Points over the external ring in the right panel of Figure 4 obey the following parametric equations: 


(A5) 


x(E,t)= A cos E + x c (t) (Bl) 

y{E,t) = BsinE + y c (t) 

where E is the eccentric anomaly, x c (t),y c (t) are the instantaneous coordinates of the planet center, and t is an 
arbitrary parameter (time for instance). 

Contact positions are those for which the distance D to origin, 


obeys two conditions: 


D 2 (E , t) = [AcosE + x c (t)] 2 + [B sin E + y c (t)] 2 , 


(B2) 


or explicitly: 


dD 2 

~dE 


= 0 

Ei,U 


D 2 (E i ,t i ) = R 2 


(B3) 


A 2 cos 2 Ei + B 2 sin 2 Ei + x 2 {ti) + y 2 (ti) + 2 [Ax c (t)cosEi + B y c {ti ) sinE] =0 (B4) 

( B 2 — A 2 ) sin Ei cos Ei — Ax c (ti) sinE + By c {ti) cos Ei — R 2 = 0. 

The solutions to these trigonometric equations, altough programmable, are not expressable in closed form. 

The “empirical” formula in Eq. (5) was obtained after the approximations B and | cos Ei\ ~ 1. This only works 
for cosir ~ 0 and 9 ~ 0. We have verified, however, that contact times estimated with Eq. (5), ( 6 ) and (7) are off by 
<1% in the case of contacts 3 and 4 for most ring projected orientations (provided 9 > 0) and <10% for contacts 1 
and 2. In the cases when 9 > 60° or cos i > 0.5 a numerical solution for Eqs. (B4) is required to attain an acceptable 
precision. 
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